Critical exponents of steady-state phase transitions in fermionic lattice models 
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We discuss reservoir induced phase transitions of lattice fermions in the non-equihbrium steady 
state (NESS) of an open system with local reservoirs. These systems may become critical in the 
sense of a diverging correlation length upon changing the reservoir coupling. We here show that 
the transition to a critical state is associated with a vanishing gap in the damping spectrum. It 
is shown that although in linear systems there can be a transition to a critical state there is no 
reservoir-induced quantum phase transition between distinct phases with non-vanishing damping 
gap. We derive the static and dynamical critical exponents corresponding to the transition to a 
critical state and show that their possible values, defining universality classes of reservoir-induced 
phase transitions are determined by the coupling range of the independent local reservoirs. If a 
reservoir couples to N neighboring lattice sites, the critical exponent can assume all fractions from 
1 to 1/{N -I). 



I. INTRODUCTION 

Experiments with cold atoms allow unmatched control 
over quantum systems giving access to a number of in- 
teresting many-body effects [H. A prominent example 
are quantum phase transitions [2] in the ground state of 
many-body Hamiltonians. On the other hand the tools of 
quantum optics allow to control the interaction with the 
environment and systems can be prepared in the non- 
equilibrium steady state (NESS) of an open dynamics 
instead [3-6]. In general the steady state of an open 
system is very different from the ground state of the cor- 
responding Hamiltonian or even from thermal states and 
thus driving a system into that state is uniquely different 
from the canonical treatment of decoherence as a per- 
turbing effect on ground states As the steady state 
or more generally the stationary space is an attractor of 
the nonunitary evolution, it is robust against further de- 
coherence. Such a scheme, where the stationary state is 
a projector, i.e. a dark state of the nonunitary evolu- 
tion, has recently been proposed e.g. by Diehl et. al. 
to prepare exotic states such as the Kitaev edge modes 
[H or BCS-type states [9]. The time evolution within a 
higher dimensional dark-space can correspond to inter- 
esting effective many-body Hamiltonians. Examples in- 
clude the recently observed dissipative Tonks- Girardeau 
gas of atoms [10] or corresponding proposals for photons 

As there is a growing theoretical and experimental in- 
terest in engineered open systems, see e.g. the recent 
review [HI, we want to discuss in the present paper the 
analogue to a quantum phase transition in Hamiltonian 
systems in the NESS of an open system, described by 
a Lindblad Liouville operator and induced by changing 
reservoir couplings. In Hamiltonian systems a quan- 
tum phase transition results from the competition of 
two non-commuting parts of a microscopic Hamiltonian 
H = Hi -\- gH2 with different symmetries upon chang- 
ing their relative strength g [2]. A transition between 
two distinct quantum phases occurs at a critical value 
g = Qc and can be identified by a non-analytic behavior 



of an order parameter. Furthermore at the critical point 
the excitation gap /S.h^ i.e. the energy difference be- 
tween the lowest excited state of H and the ground state, 
closes. At the same time certain correlations become in- 
finitely long-ranged indicated by a diverging correlation 
length ^. Criticality induced by reservoirs in open many- 
body systems has recently been discussed by Eisert and 
Prosen [12] for free systems, where criticality was defined 
in terms of a diverging correlation length. Here we will 
reexamine free fermionic lattice models as generic models 
to study noise-induced quantum phase transitions. As 
analogue to the excitation gap /S.h in unitary systems 
we will consider the damping gap A of the Liouvillian as 
indicator of a critical point for reservoir-induced phase 
transitions. We will show that if the manifold of reser- 
voir parameters is of dimension d the values for which 
the damping gap closes and the system becomes critical 
is generically of dimension d — 2 for free systems. As 
a consequence the gapped phases are always connected 
in parameter space and there are only isolated critical 
points or manifolds for reservoir-induced quantum phase 
transitions of linear systems. 

An important aspect of quantum phase transitions in 
unitary systems is that close to criticality, microscopic 
details of the interaction become irrelevant leading to 
universal behavior and allowing to classify phase tran- 
sitions according to universality classes. These classes 
are defined by common critical exponents A and dynam- 
ical critical exponents characterizing the divergence of 
the correlation length and the closure of the excitation 
gap when approaching the critical point: 



We here argue that similar universal exponents can be 
derived for reservoir-induced transitions where the possi- 
ble critical exponents depend solely on the coupling range 
of the independent local reservoirs. 

After introducing linear fermionic lattice models cou- 
pled to local reservoirs and summarizing the general 
methods to derive the NESS of these systems in section 
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II, we discuss a specific realization of such a model with 
reservoirs coupling to two neighboring lattice sites in cold 
atom systems in section III. Using this example as illus- 
tration we then formulate the key questions addressed 
in the remainder of this paper: Under what conditions 
do reservoir induced phase transitions occur and to what 
universality classes can they be attributed to? We will 
then show in section IV for general reservoir couplings 
that noise-induced criticality defined by a diverging cor- 
relation length occurs at points in parameter space where 
the damping gap, i.e. the spectral gap of relaxation rates 
of the Liouvillian dynamics closes. In contrast to Hamil- 
tonian systems, where the parameters that drive a phase 
transition are real, in Liouvillian system these parameters 
are complex. We then show that for linear models crit- 
icality occurs only at isolated singularities or parameter 
manifolds of dimension d — 2, where d is the dimension of 
the reservoir parameters. An immediate consequence of 
this is that there is only one connected phase with a non- 
vanishing damping gap. We derive the possible critical 
exponents, defining universality classes of noise-induced 
phase transitions and show that they are simple frac- 
tions determined by the number of sites that couple to 
the same reservoir. Finally we establish a simple relation 
between the critical and dynamical critical exponents. 



II. FREE LATTICE FERMIONS WITH LINEAR 
RESERVOIR COUPLING 



As a generic system we study a free and translationally 
invariant fermionic chain described by the system Hamil- 
tonian H. In addition to the unitary evolution there is a 
non-unitary interaction with reservoirs which is assumed 
to be Markovian. This allows for an effective description 
of the dynamical equation for the density operator g via 
Lindblad operators 



d 



operators 

2L 

= ^ {H)j^kWjWk, (3) 

L 

= '^{l)f^,j,l^2j-l + (0/i,J,2^2j. (4) 

J = l 

In the following we consider local reservoirs only. I.e. we 
restrict the reservoir coupling to a finite number of N 
adjacent sites. To explicitly take into account transla- 
tion invariance we reformulate the Hamiltonian and the 
Lindblad generators 



j \m,n / 

N-1 

m=0 



(5) 
(6) 



where we have introduced the operator r^, which shifts 
a local operator by j lattice sites. One recognizes that 
the Lindblad generators contain 2N independent com- 
plex parameters Sm = fm^^^"^ and qm = z>rne*^"^. Only 
27V — 1 are relevant however in the context of reservoir in- 
duced phase transitions as one coefficient can be pulled 
out and only determines the overall time scale of the 
damping. Thus we may set without loss of generality 
50 1. 

The dynamical equation ([1]) contains only quadratic 
terms in fermionic operators, thus the NESS is Gaus- 
sian [l^ and is fully described by the correlation matrix 



-Ti{gQ(wjWk 



WkWj 



(7) 



\ (^^M^^M ~ O}) ' *^ ^ ^exp(^ - ^OjTjkOk^ . The dynamics generated by Eq.© 



Higher order correlation functions can be calculated us- 
ing Wick's theorem [1^. Tjk is antisymmetric, has 
purely imaginary entries and is directly related to the 
Grassman representation of the steady state w{go,6) = 



In the following we will use Major ana operators 



take the simple form of a linear matrix differential equa- 
tion for the correlation matrix F Il2l 



W2j-1 



^3 •) 



W2j = i{Cj 



4)' 



(2) 



dt 



X' F + FX - r. 



(8) 



instead of the more familiar fermionic creation and an- 
nihilation operators Cj and Cj, as the Major ana opera- 
tors allow for an easy representation of the steady state 
of the considered linear systems. They are analogues 
of the bosonic position and momentum operators, are 
hermitian and fulfill a simple anti-commutation relation 
{wj^Wk} = '^Sj^k- The signature of a free system is the 
bilinearity of the Hamiltonian and the linearity of Lind- 
blad generators when expressed in terms of Majorana 



The matrices X and Y represent the Hamiltonian and 
the Lindblad generators X = -AiH - {R^ R*) and Y = 
2i{R — i?*) with the reservoir matrix R = ^ Ifi^l*^- 

We are only interested in the steady state of the system 
given by the solution of the Lyapunov Sylvester equation 



X^Fo + FoX = r. 



(9) 



For translational invariant and infinitely large systems, 
the correlation matrix of the steady state is circulant 
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and can be represented by its Fourier transform, called 
the symbol function 7((/^), which is a 2 x 2 matrix corre- 
sponding to the two different types of Major ana species 
with even or odd indices: 

li^) = ( ^''\^\ ) (10) 

^^^^ \l2l{(p) 722 (V^) J ^ ^ 

where 

(11) 

and analogously for other types of Majorana operators. 

In terms of the 2x2 symbol functions the Lyapunov 
Sylvester equation reads 

^(-^)^7(^) + 7(^)^(^) = y{^)- (12) 

The matrices x{ip) and y{(^) are calculated in correspon- 
dence to the X and Y matrices of the finite size model. 




FIG. 1. (Color online) Fermions with four internal states 
are trapped in state dependent optical lattices in the above 
configuration to realize the reservoir coupling of two sites in 
the 1^) lattice. The shallow potential for atoms in internal 
state |r) acts as the reservoir, which is coupled in a controlled 
way to 1^). The effective pumping is driven by a laser on the 
|ei),|r) transition, whereas the decay is driven on the |e2),|^) 
transition. 



III. A QUANTUM-OPTICAL EXAMPLE: 
LATTICE FERMIONS WITH TWO-SITE 
RESERVOIR COUPLING 

In order to illustrate the physical context of the present 
paper let us consider fermionic atoms with four internal 
states 1^'), |r), |ei), |e2) in state selective, optical lattice 
traps as shown in Fig. [TJ The lattice for atoms in states 
\g)-, 1 62) is given by an optical standing wave and has 
lattice constant ^. The ground state lattice is shifted 
by half a lattice constant compared to the other two. 
Atoms in internal state |r) feel a very shallow poten- 
tial and are delocalized compared to the tightly confined 
atoms in other internal states. The transition |r) ^ |ei) 
is driven by a laser field with Rabi frequency (^i, whereas 
1^) ^ 1 62) is coupled to a laser with Rabi frequency 1^2- 
Spontaneous decay occurs from the two excited levels into 
the metastable ground states \r)^\g). The optical lattices 
are deep so that we can assume that only adjacent Wan- 
nier wavefunctions (/pf(x),/i G 61,62 and ip^ or ^^j^i are 
optically coupled by the laser fields. Using this setup 
non trivial pump and loss processes into and from the 
metastable states \g) are realized. Atoms in the shallow 
|r) potential act as a reservoir for the optical transitions. 

First of all atoms from the reservoir |r) are pumped 
via 1 61) into a superposition of atoms in neighboring lat- 
tice sites with fixed relative phase. To see this we note 
that the spontaneous emission from \ei) to \g) can be 
described by the interaction Hamiltonian 

n^nt = j dx'^(^g^a^¥g{x)i!,,{x)e'^^^ (13) 

k 

where is the annihilation operator of the electro- 
magnetic mode with wavevector /c, and g^ is the corre- 
sponding coupling matrix element. Using the decompo- 
sition of the fermionic fields in internal states \g) and 



\ei) into the Wannier basis ^ g{x) = ^j^^j{x)cj and 
^ei(^) = E,^f (^)el,j yields 

Uint = ^Y1 [Sk^k^j + ^j2^4+l) + ^'^] ' 

j k 

k ^^^^ 

where 77^^^^ = J dx cp^^ {x)(p^{x)e'^^''^ and 77^2^ = 
J dx (f^j^ {x)(p^j_^-^{x)e'^^''^ denote the Frank-Condon fac- 
tors corresponding to the transitions j ^ j, j + 1. Due to 
the exponentially decreasing Frank- Condon overlaps all 
transitions with / 7^ j, j + 1 can safely be neglected. As 
the Wannier functions in a deep optical lattice are well lo- 
calized the products cp^^ {x)(p^{x) and cp^^ {x)(p^_^-^{x) are 
well localized functions at positions Xj±a/4 with a ~ A/2 
being the lattice constant. Thus (fT4|) can be rewritten as 

Hint = E E [dkH^^^i (4 + 4+1^''^"/') + h.a: , 
i k 

(15) 

where the (real) parameter v can be tuned by shifting 
the position of the \ei) lattice, relative to that of \g). 
Coupling of the many motional states in |r) with a laser 
to 1 61) leads after elimination of the vacuum modes to an 
optical pumping that can be described by independent 
Lindblad generators 

i™ = x(4 + ^+i) (16) 

with X ^7 which describe the coupling to two adjacent 

lattice sites. Note that the relative phase term e*^^^/^ in 
Eq. (p!5|) vanishes after averaging over the vacuum modes 
up to the first order in kQa/2. 

We now show that optical pumping from \g) via I62) 
leads to a loss of fermions in all superpositions of neigh- 
boring sites except for one dark mode. The correspond- 
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ing Hamiltonian describing the coherent part of the in- 
teraction reads 

^int,2 = ^ iVjiCj + fjj2Cjj^i) + h.a. (17) 



J dx (p^^ {x)(pj{x)e'^^''^ and fjj2^ 



where 7)jJ^ 

/ dx (p^j^ {x)(p^j_^-^{x)e'^^''^ where Qx = q - and q is the 
wave vector of the laser corresponding to 0.2- Note that 
since the wave vector of 1^2 is weh defined fjji and fjj2 
differ in both amphtude and phase. Considering a fast 
subsequent decay from |e2) finally gives rise to an opti- 
cal pumping out of state \g) described by independent 
Lindblad generators 



(18) 



Without loss of generality we can set 7 = 1 which fixes 
the overall time scale of the process. Then the free pa- 
rameters of the Liouvillian are the amplitudes x u 
and the phase g. 

Solving the Lyapunov- Sylvester equation (|9|) we find 
for the symbol function of the correlation matrix 



(19) 



d{(p) \-ni2{(p) nii{(p) J ' 
where 

^11 (^) = 4zi^x^ (1 + + 2i/ cos (p) sin g sin (/?, (20) 
^i2(^) = (1 + + cos^ cos (/:))^ (21) 

00 A. / \2 

— sin ^sin (p — x (1 + ^ -\-2ucos(p) 

and 

d{ip) = [(1 + i^^) + +2^(x^+cos^)cos(^]^ 

- 4i/^ sin^^sin^ (p. (22) 

The Fourier-transform according to Eq. (pT]) yields the 
correlations of Major ana fermions. From the symmetry 
of j{(p) it is immediately clear that as expected only nor- 
mal correlations of fermionic creation and annihilation 
operators are non zero. Following Ref.[12] we can define 
criticality by a diverging correlation length 



lim 



log 



'SS 



(23) 



The correlations (cjcj) become infinitely long ranged for 
u = Uc = 1 and g = gc = and any non- vanishing value 
of X- In Fig- El^a) the dependence of ^ on the phase g — gc 
is shown for u = 1. One recognizes a linear dependence 
and the same holds for the dependence on u — Hence 
the critical exponent for this example is A = 1 . The same 
critical exponent has been found by Eisert and Prosen in 

The NESS with correlation matrix Fq is an attractor 
of the dynamics, i.e. small deviations ^F from it will 




(b) 



1.6 



A 0.8 




FIG. 2. (Color online) For the quantum- optical example the 
inverse correlation length shown in (a) shows a linear behav- 
ior around a critical point at = 0, z^c = 1- Near the critical 
point the damping gap A, shown in (b), approaches zero with 
exponent 2. Both length and time scales diverge in the vicin- 
ity of the critical point. 



decay to zero after some time. An important quantity 
is the smallest decay rate for such deviations as it de- 
fines the time-scale of decay back to the stationary state. 
Since close to the critical point correlations become in- 
finitely long ranged one expects that the time scale for re- 
establishing long-range order after a small perturbation 
tends to infinity. This corresponds to a closure of the 
damping gap A, i.e. the smallest non-vanishing eigen- 
value of the real part of the Liouvillian. The damping 
gap A is here the direct counterpart to the excitation 
gap Ah in unitary systems. Staying in the manifold of 
Gaussian states, the dynamical equation for 6r reads 



—ST = X^ ST^STX. 
dt 



(24) 



Thus one has to consider the eigenvalues of X, i.e. 
R -\- . One again finds that these eigenvalues vanish 
for arbitrary values of x for = = 1 and g = gc = ^^ 
Fig. Efb) shows the damping gap for the present ex- 
ample. One recognizes a closure at the isolated point 
Sc = 1^0^^^'' = 1 with a quadratic dependence in the im- 
mediate vicinity. Thus the dynamical critical exponent 
is z = 2. 

In this particular example we have seen that a critical 
point is associated with both a vanishing of the damping 
gap A and a diverging correlation length. One recog- 
nizes furthermore that although the parameter space of 
the reservoir coupling is two-dimensional (disregarding 
the irrelevant parameter x), there is only a single point 
where the system is critical. The parameter region with 
nonvanishing A is a connected manifold and thus there 
is only one distinct gapped phase. Continuous changes 
of the parameter s = ve^^ in the complex plane circum- 
venting the critical point Sc = 1 will smoothly connect 
the entire gapped region. 

The present example gives rise to a number of ques- 
tions: Is a diverging correlation length always connected 
to a vanishing damping gap? What are the possible uni- 
versality classes, characterized by the critical exponents 
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A and the dynamical critical exponents zl In the follow- 
ing we will discuss these issues for linear fermion models 
with general local reservoir couplings. 



IV. FREE LATTICE FERMIONS WITH 
GENERAL LOCAL RESERVOIRS 



Let us now consider fermionic lattice models with 
reservoirs that couple simultaneously to N adjacent lat- 
tice sites. A schematic representation of our model is 
shown in Fig. [3l To simplify the calculations and the final 
expressions we restrict ourselves to Lindblad generators 
which only contain a single type of Major ana fermions. 
A generalization is however straightforward and all con- 
clusions hold. In particular we consider only Majorana 
operators of the first kind {w2j-i) in the generators as 
they are invariant under the exchange of creation and 
annihilation operators. The dissipative dynamics is de- 
coupled from the even Majorana modes, leading to de- 
generacies in the NESS. This degeneracy can be lifted by 
a free, translation invariant Hamiltonian and the addi- 
tion of this unitary term does not change the properties 
of the uneven Majorana modes. In this case the steady 
state equation ([12]) has a trivial solution in terms of the 
symbol function r ((/:?), which represents the reservoir cou- 
pling 



J^mJ^ne-*'^'™-")e*(^''"-S"\ (25) 



r(y) - r(-(/;) ^ 
r((/?) + r(-y)) 



(26) 



Apart from the fact that the Hamiltonian guarantees 
its uniqueness, the NESS is independent of the Hamil- 
tonian details and therefore possible critical features of 
the ground state of 1-i are not recovered in the steady 
state. Moreover ([26]) does in general not correspond to a 
pure state. The particle-hole symmetry of the Lindblad 
generators leads to a mean occupation of \ in the NESS. 
The eigenvalues of the circulant correlation matrix F are 
given by the entries of the symbol function 7((/:^) which 
is positive. They are bounded between [—1,1], which is 
in contrast to pure states for which it can be shown that 
all eigenvalues must be ±1. 

By inverse Fourier transform we calculate real-space 
correlations in the steady state 



1 r^'' 



(27) 



The integration can not be carried out in general but 
we can understand the characteristic properties for large 
spatial distances d by using arguments of complex calcu- 
lus. To this end we rewrite the symbol 7((/?) as a function 




FIG. 3. (Color online) (a) Schematic representation of a linear 
chain translationally invariant coupled to independent reser- 
voirs, (b) The single reservoir is characterized by 2N real 
parameters as indicated for a three site reservoir coupling. 



in the complex plane using the substitution e*' 



711 (^) 



{wiwi^2d) = ^ Kesa [z"^ ^lii{z)] , 

aeSi 



(28) 
(29) 



where Res^ denotes the residues inside the unit circle , 
{\z\ < 1). The residue is non-zero only in singular points 
of ^{z). Because numerator n{z) and denominator d{z) 
of (|29]) are holomorphic, only zeros of the denominator 
inside the unit circle contribute to the correlations. The 
symbol function has at most simple poles and it has been 
pointed out in [l2| , that the zero closest to the unit circle 
is relevant for the large d behavior. As the denominator 
d{z) of Eq. ([26|) is just the symbol function of the matrix 
X, zeros of d{z) on the unit circle correspond to both a 
diverging correlation length and a closure of the damping 
gap A, making both definitions of noise-induced critical- 
ity for linear models identical. 



A. conditions for criticality 

The critical points, i.e. the singularities of the sym- 
bol function are the roots of d{z) on the unit circle. 
For general reservoir couplings to multiple sites, ex- 
plicit expressions for the roots of d{z) are hard to ob- 
tain and we need a different criterion to find the criti- 
cal parameters. On the unit circle the denominator d{z) 
is strictly non-negative because it can be rewritten as 

configuration {sj} = {z^je*^^ } of the complex parameters 
of the Liouvillian leads to a critical behavior if a zq on the 
unit circle exists, such that the individual sums inside the 
absolute value vanish for zq and its complex conjugate. 
This gives a pair of implicit equations 



AT-l 

1 + ^0 - 



N-1 

E 

i=i 



J* _ 



(30) 
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where we have used that without loss of generahty sq can 
be set equal to unity. Apparently reservoir couplings to 
a single site (A^ = 1) cannot induce criticality, however 
couplings to N > 1 sites may. For a given zq, 2{N — 2) 
out of the 2{N — 1) real parameters I'j^gj can be chosen 
arbitrarily. As there are only a finite number of roots zq 
the non- trivial complex solutions {sj = VjC^^^} to these 
equations are limited to a (i — 2 dimensional manifold in 
the d = 2 (A' — 2) dimensional parameter space. As a con- 
sequence there can never be two extended, non connected 
regions in parameter space with a finite damping gap A. 
Thus linear systems can become critical, but there are 
no reservoir-induced phase transitions between distinct 
gapped phases. 



B. correlation length and critical exponents 

In the vicinity of the critical point the behavior of ^ is 
determined by the leading-order exponent A of the sin- 
gularity, which itself is determined by the properties of 
the denominator d{z) in Eq. (|28l) . If a zero of d{z) ap- 
proaches the unit circle from the inside, ^ diverges which 
corresponds to a phase transition to criticality. Let zq be 
the closest singularity to the unit circle, then the corre- 
lation length is given by 



-In 1^0 1 



1 



l^ol- 



(31) 



In the following we will analyze the dependence of ^ on 
the system parameters in the vicinity of such singularities 
and determine the corresponding critical exponents. 

For two-site coupling the implicit equations can have a 
nontrivial solution and we find that the generator Ltwo = 
uqc^^^wi -{-uie'^^^ws leads to a critical NESS for go— gi G 
ttWj and = Analyzing the behavior in the vicinity 
of the critical points we find a critical exponent of A = 1 
in agreement with the results of [l2^ . The question arises 
if this value is the only possible one in free fermionic 
lattice models. To answer this question let us consider 
a more general reservoir with coupling to three adjacent 
sites: Lthree = ^two + ^^2e'^2^^;5. 

A possible solution of the implicit equations (|3Q|) is 
= ui = V2 and g^ = = 0,^2 = Under 

variation of for example gi the critical exponent is here 
again A = 1. However, there is another solution go = 
gi = g2 and J^o = ^^1 = ^2' In this case the critical 
exponent under variation of gi is different and given by 
A=i. 

In Fig. m we have plotted the two-point correlation 
functions for sites separated by a distance d in depen- 
dence of the phase g that drives the transition for the 
case of a three-site reservoir coupling. The left part of 
the figure corresponds to the A = | case, whereas the 
right part corresponds to a transition with A = 1. The 
analytic behavior of the inverse correlation length around 
the critical points, shown in the lower parts of the plots, 
clearly distinguishes the two cases. 



(a) -0.5 
-0.25 




0.25 
0.5 



n 



15 30 45 60 15 30 45 60 
(b) ^ (d) ^ 



-2 
-4 
-6 
-8 
-10 



0.4 



-0.4 0.4 

9 -9c 




-0.4 0.4 

9 -9c 



FIG. 4. (Color online) (a) and (c) show two site correla- 
tions \wiwi-^2d\ on a logarithmic color scale in dependence of 
the spatial distance d and a parameter change around critical 
phase Qc. The corresponding inverse correlation length is 
plotted in (b) and (d) and reveals the critical exponent. The 
left side is the A = | transition and the right side the A = 1 
transition for the chain with three site reservoir coupling dis- 
cussed in the text. 



One recognizes that the amplitude of the component 
with long-range correlations vanishes as one approaches 
the critical points. More precisely for the case of the 
smaller critical exponent, A = 1/2, all non-local corre- 
lations vanish, while for the case of the larger critical 
exponent, A = 1, some component with finite correlation 
length remains also at the critical point. 

In the following we want to get some general in- 
sight into what are the possible critical exponents of a 
reservoir-driven phase transition. As we are interested in 
the behavior in the vicinity of the critical point, we need 
to expand the denominator d{z) of the symbol function in 
terms of the relevant system parameter {^j}, {i^j} around 
their critical values. Since we do not have an explicit ex- 
pression for the roots zq of d{z), we need to do this in an 
implicit way, i.e. expanding d{z) both in terms of z — zq 
and in their explicit dependence on {^j}, {^^j}. 

Due to the positivity of d{z) on the unit circle, all 
first order partial derivatives with respect to z as well 
as to the system parameter gi^ and Uk must be zero at 
the critical point. To find the leading order expansion in 
z — Zq we thus evaluate higher order partial derivatives 
with respect to z using the implicit equations (|3Q|) 



d^d 



0^ = ^ 1^ 



K93-91) 



[j - - I - l)z- 



3-1 



3,1 



+{i-m-3-i)z'-' 



(32) 



If the second order derivative is nonzero at z = zq, we 
can stop at this level. On the other hand, the second 
order derivative vanishes if a second pair of independent 
implicit equations is fulfilled, which can easily be read 
off from (|32]) . This procedure can be continued and each 
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term 



equations 



which is zero, yields a new pair of imphcit 



(33) 



Here we have used that the first non- vanishing derivative 
must be an even one. Let us assume that ah derivatives 
in z vanish up to order 2M — 1. Then it can be shown, 

that mixed derivatives of the type d(^^^^g^^d'^d vanish 

for ah m < M. Thus what remains are the second order 
partial derivatives with respect to Vk or gk- Second or- 
der partial derivatives in the same parameter are always 
non zero on the unit circle (except for trivial cases), as 



dUiz) 



2- 



2i/|. Thus we can write 



the power expansion of d{z) in the following general way 



x) ^ C2X' 



-2 ' -C,^mZ^ + 0{z^+^) 



(34) 



where C2 and Cq ^^^^ constants. Here x is 

a linear combination of parameter variations from the 
critical values h = 9k - Qkc and i>k = J^k - i^kc, and 
z = z — Zq. The lowest non- vanishing contributions de- 
termine the critical exponent and therefore M is the min- 
imal M of all parameters included in x. The zeros zq are 
algebraic functions of the system parameters and we can 
therefore write z ~ + O^x'^^^). At least two terms in 
the expansion must be of the lowest order and therefore 
we find A = along the line x if the first M implicit 
equations are fulfilled. We see that all possible critical ex- 
ponents are the inverse of integer numbers. The smallest 
possible critical exponent is determined by the maximum 
M, which is just given by M. 

We now can relate the M to the number N of adja- 
cent sites coupled by each local reservoir. It is clear that 
equations (|33|) are linearly independent for different m. 
On the other hand only N equations can be independent 
for a finite reservoir coupling. This proves that M can 
be at most — 1 or all orders vanish, in which case the 
symbol function must be zero everywhere and the system 
is not critical. We conclude that if the reservoir coupling 
is restricted to N sites, the critical exponent is out of a 
bounded set of fractional numbers 



AG U 



1 1 
''2' 3'"' 'TV 



(35) 



This is the main result of the present paper. 

The corresponding Taylor expansion of the numerator 
n{z) in a critical point is of higher order than the de- 
nominator. Therefore the amplitude of the critical corre- 
lations, vanishes as ^ diverges. This is seen in the graphs 
of Fig. m The remaining non-local correlations, visible 
for example in part b) of that figure are due to additional 



singularities inside the unit circle. Only in the case of the 
minimal critical exponent A = 7^^, these singularities 
cannot exist due to fundamental laws of algebra. In this 
case the NESS is completely mixed in the critical point. 

Another result that can be drawn from our analysis is 
the dimensionality of the critical parameter space. Criti- 
cal points have to fulfil the set of equations ([33|) and for a 
given critical exponent the corresponding dimensionality 
is given by 



dim(PA) = 2 A- - 1 



(36) 



It is clear that the critical points are always a zero mea- 
sure subset of parameter space, but they are not neces- 
sarily isolated points, with the exception of critical points 
with minimal exponent for the given configuration, which 
are always singular. 



C. spectral gap of relaxation rates and dynamical 
critical exponent 

For unitary lattice models it is well established that 
the presence of a finite gap in the excitation spectrum 
leads to a finite correlation length, while the transition 
to criticality is associated with a vanishing gap [15]. In 
the following we want to establish a corresponding rela- 
tion for reservoir driven phase transitions and discuss the 
dynamical critical exponents. 

The relaxation rates of the system are determined by 
the homogeneous part of (j8j) and therefore the damp- 
ing matrix X. More precisely the damping matrix X 
describes the dynamics on the submanifold of Gaussian 
states, whereas the full system is spanned by the Liou- 
ville Operator C in equation ([1]). The trace preservation 
of the Lindblad dynamics requires C to have at least one 
eigenvalue with vanishing real part. If this zero eigen- 
value is unique, the gap in the real spectrum sets the 
slowest relaxation rate for arbitrary initial states. The 
gap A of X thus gives an upper bound for the gap of £, 
and the gap in the full damping spectrum must vanish 
as one approaches the critical points. The eigenvalues 
of X are purely real, when neglecting the Hamiltonian 
contributions and strictly negative. In the translation 
invariant system the eigenvalues are given by the symbol 
function —r{(f) — r(— (/p), which we have identified before 
as relevant for the correlation length. In the vicinity of 
a critical point, the slowest relaxation rate, defining the 
spectral gap of relaxation, is determined by the roots Zj 
of (i(2:, {^j}, {^j}) closest to the unit circle 



id{z,{gj},{vj}). 



(37) 



A mmi^i^ic 

Therefore the dynamical exponent is immediately related 
to the critical exponent A. The exponent of the diver- 
gence however must be modified by the number of roots 
He-, that merge at the same point on the unit circle when 
the Liouville parameters approach their critical values: 



\9-9c 



(38) 
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The number of roots i^c is thus identical to the dynamical 
critical exponent z. For the minimum critical exponent 
A = 1/{N — 1) all 2(A^ — 1) complex roots merge simul- 
taneously on the same point and thus 

zX = 2, for A = Amin (39) 

Moreover for all examples we have considered we found 
that the damping gap closes as a quadratic function, 
which suggests that (|39|) is more general. 

V. SUMMARY 

To summarize we have analyzed the non equilib- 
rium steady state of translation invariant chains of free 
fermions coupled to local Markovian reservoirs described 
by linear Lindblad generators. Such couplings can be 
generated e.g. in ultra-cold atomic lattice gases as we 
have shown for the example of the two site coupling. 
A general expression for the correlation matrix of the 
NESS can be obtained using a symbol function ansatz. 
We showed that under certain conditions the NESS goes 



into a critical state upon changing reservoir parameters. 
The critical state is characterized by a simultaneous di- 
vergence of a correlation length and a critical slow-down 
of relaxation i.e. a closing of the gap in the damping 
spectrum. We showed that the dimension of the crit- 
ical parameter space is at most d — 2, where d is the 
dimension of the full space of reservoir parameters. As 
a consequence all gapped phases are smoothly connected 
and although there is a transition into a critical phase 
there is no reservoir-induced quantum phase transition 
between distinct gapped phases for linear models. The 
transitions to a critical state can be classified by the lead- 
ing order exponent in the dependence of the inverse cor- 
relation length on the system parameter that drives the 
transition. We have shown that this critical exponent 
must be the inverse of an integer between 1 and — 1, 
where N is the number of sites coupled by a single reser- 
voir. Furthermore a general expression for the dynamical 
critical exponent, describing the closure of the damping 
gap was derived. 
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